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An intensive investigation inlo the cause of anomalous behavior commonly observed in impedance eduction 
techniques is performed. The investigation consists of grid refinement studies, detailed evaluation of results at 
and near anti-resonance frequencies, comparisons of different model results with synthesized and measured 
data, assessment of optimization techniques, and evaluation (if boundary condition effects. Results show that 
the root cause of the anomalous behavior is the sensitivity of the educed impedance to small errors in the mea- 
sured terminal ion resistance at frequencies near anti-resonance or cul-on of a higher-order mode. Evidence 
is presented lo show (ha I the common usage of an anechoic, plane wave termination boundary condition in 
ducts where the “true” termination is reflective may act as a trigger Tor these anomalies. Replacing the exit 
impedance boundary condition by an exit pressure condition is shown to reduce the anomalous results. 
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iinite element system matrix, vector of source effects 

complex mode coefficients, complex constant for defining transcendental equation 

ambient sound speed, mean static temperature 

material derivative, 2D gradient operator, 2D Laplace operator 

depth of test liner, axial location of duct termination, height of duct (0.051 meters) 

transcendental equation, average (uniform) Bow Mach number 

excitation frequency, angular frequency 

index of summation, free space wavenumber 

unit imaginary number, number of wall measurement points 

leading edge of test liner, trailing edge of test liner 

axial propagation constant, mode eigenvalue 

mean static pressure in duet, mciin static density in duel 

number of points in axial direction, number of points in transverse direction 

acoustic pressure field, ax i til particle velocity, transverse particle velocity 

acoustic particle velocity vector, mean flow velocity vector 

transverse, spanwise, and axial coordinate 

time, location of a lower wall pressure measurement 

ratio of specific heats, normalized acoustic impedance 

normalized resistance, normalized reactance 

vector of unknown node coefficients, wall objective function 

absolute value of complex quantity, vector dot product 
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Subscripts: 
0, wall 
ave, ini lial 
s,exit 
FEM . M 


mean flow quantity, wall quantity 

spat ially av eraged q rum l i l y , initial as si goed v al uc 

source plane quantity, exit plane quantity 

quantity computed from finite element analysis, quantity measured on lower wall 


I, Introduction 

R ecently, the NASA Langley Liner Physics Group conducted a detail led assessment of a number of propagation 
codes, one of which is used as the basis for the InRC impedance eduction technique. This detailed assessment 
shed valuable light on potential causes for previously documented anomalous behavior; i.e., educed impedances that 
did not match the expected trends. This anomalous behavior has been observed with the single mode model, 1 the 
convecied Helmholtz model, 2 and the linearized Euler model, 5 and therefore does not appear to be caused by the 
choice of model. These anomalies have been observed at frequencies near anti -resonance and or cut-on of higher- 
order modes. Although previous studies 5 have documented potential explanations for the anomalous behavior, the 
root causes have not been fully investigated. Hie purpose of this paper is to launch an intensive investigation into the 
cause of the anomalous behavior. To this end, data are evaluated across wide frequency and Mach number ranges, such 
that detailed trend studies can be observed at frequencies near anti -resonance and/or cut-on of higher-order modes. 

II. Statement of Problem 

T ins analysis considers acoustic wave motion through a uniform flow, constant area, rectangular duct as shown 
in Fig. 1. All walls of the duct (with the exception of the portion containing the liner) arc considered rigid mid 
impervious to sound. The lining is locally leading and the unknown normalized impedance of the portion of the 
lop w;ill containing the lining is denoted by ^ wa ii- Throughout this report, all impedances are normalized with the 
characteristic impedance, poCo, of the air flowing in the duet. A right-handed Cartesian coordinate system (x, y,z) 
is used, with the z-axis pointing dong llie axis of the duct, the y-axis perpendicular to the duel sidewalls, mid the 
x-axis perpendicular to the top and bottom walls. Consistent w ith most experimental apparatuses and aircraft engine 
nacelles, the source and termination planes are located in rigid wall sections of duct. The excitation frequency is 
assumed below cut -on of higher-order modes, so that the span wise direction (i.e., the y direction) supports only plane 
waves. Additionally, it is assumed that the acoustic pressure on the wall opposite the test liner, Pm(z ) 7 is measured at 
a number of axial stations, z = z/, and that the duct termination is near mi echoic. 

The problem al hand is to determine the unknow n impedance of the test liner given a measurement of the source 
plane conditions, the acoustic pressure on the bottom wall opposite the lest liner, Lhc termination plane impedance, 
and the mean flow held. More specifically, it is desirable to thoroughly evaluate I^aRC's two impedance eduction 
methodologies [the convecied Helmholtz model (CHF,) 2 and the linearized Euler model (LEE) 4 ] on the same liner, 
and over a range of frequencies and flow Mach numbers. The goal of this evaluation is to determine the source of the 
impedance eduction anomalies discussed in the previous section. This evaluation is based on a conventional, perforate- 
over-honeycomb liner that is typical of that currently employed in high-bypass aircraft engine nacelles. Data used to 
perform this investigation are either synthesized from mode theory, or measured in the NASA Langley Research Center 
Grazing Incidence l ube (GIT). 


III. The CHE Impedance Eduction Model 


r j^HE 2D CHE model is based on the solution to the 2D w ave equation in a uniform flow' 


I D 2 p 


V 2 ?, 



4* ^ 


( 1 ) 


w here a lime dependence, has been assumed. In the current subsonic flow problem, the convccted wave equation 

that governs the acoustic pressure field, i.e,, Eq. (1), requires one boundary condition at the source boundary (z ~ 0), 
another at the duel terminal ion (z = L), and another at a rigid or acoustically lined wall. The source plane acoustic 
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pressure is used as the source plane boundary condition 


P = Ps (2) 

whereas, aL Lhe near anechoic duct termination, z— L , the boundary condition is written in the form 2 

$£ = - ikp (3) 

dZ T" 3<fave 

When l he termination is anechoic and located in a hard wall duct where only planar waves me cut -on, the measured 
value of ^exit is unity (re., Ccxit = I ) ■ However, when the termination in near anechoic mid located in a hard wail section 
(e t g., as considered here) the termination impedance will vary slightly from unity. The GIT was designed to be as 
close to anechoic as possible, so that Eq. (3) is expected to represent a good approximation* Measurements of the 
exit impedance in the GIT, confirm that il is near anechoic, especially at the lower values of Mach number. It should 
be noted that Eq. (3) is also valid when a nonreflecting termination is located in a soft wall duct. However, when the 
termination is located in a softwall duet, the exit impedance (i.e., may vary significantly from a value of unity. 
Several examples will be presented in the result section for which the termination is located in a softwall section of 
duct. However in theses examples, t e xii is synthesized from a mode solution. Finally, w hen the termination is highly 
reflecting but located in a hard wall duct where only plane waves are cut -on, Hq. (3) will not be valid but can be replaced 
by the more general condition 

dp ik 
dZ 

In addition to Hq, (3), the use of the termination pressure 

P = Pent (5) 

as a boundary condition at the duct termination is also investigated. The tliird and final boundary condition required 
for the CHE model is the wall imped^inec boundary condition. When written in terms of the acoustic pressure, the 
wall impedance boundary condition is 5 
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It is noted that Hq. (6) is also used along the rigid wall portion of the duct wall. However, along a rigid portion of the 
duct wall, t^ai! is set to infinity (£ wa u — °°). 


IV, The LEE Impedance Eduction Model 

I N the 2D LEE model, the governing differential equations are the linearized equations governing conservation of 
mass 

C ®15t (?) +^0? *5& + pV + V) ^ =0 (7) 

and momentum 

ypo~^+ypo (V* Vj T ^+P~^ t - + c o^P=° ( g ) 

Here, the assumption is made that the acoustic process lakes place homen tropically in an ideal gas mid that Lhe fluid is 
inviseid. Note that Hq. (7) is a scalar equation and Eq (8) is a vector equation with two components in 2D. The above 
acoustic system requires two boundary conditions at the source plane, one boundary condition at the duet termination, 
mid a third boundary condition at mi acoustically treated or rigid duel w all. 6 At the source plane, the acoustic pressure 
and the transverse particle velocity components, are specified 

W={s) 

At t lie duct termination, the boundary condition is specified in the form 


Po^oCexil 


(9) 

( 10 ) 
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and the wall impedance boundary condition is of the form 


prv = 



M aV e d\ 

ik dz) 


P 

CwallJ 


(ID 


Equation (11) is used at both a rigid and acoustically treated wall. However, at a rigid wall, the wall impedance is set 
to infinity (i.e., ^ wa ii — rc) r 


V. Numerical Solution for the Acoustic Field 

T HE numerical method chosen to obtain the solution to the acoustic field in the CHE and LEE models is the finite 
element method. Details of the implementation of the methodology are described in several papers 2 4 - 7 and are not 
repeated herein. It is noted only that there are N and M evenly spaced nodes in the axial £ind transverse directions of 
the duct, respectively. Cubic Hermile polynomial basis functions tire used to approximate the solution to the acoustic 
pressure held in the CHE model mid a weak formulation is employed to incorporate the impedance mid near mice hoi c 
boundmy conditions. On the other hand, the LEE model uses linear elements to represent the acoustic pressure mid 
particle velocity fields within each element, mid the wall impedance mid exit impedance boundary conditions are 
satisfied by constraining the nodal degrees of freedom. For both the CUE and LEE models, the finite element method 
results in a discrete system of equations of the form 


[A]{<H} = {F} (12) 

where [A] is a sparse, complex, asymmetric matrix whose order is 4 MN mid 3 M/V, respectively, for the CHE mid LEE 
models. Further, the vector {d>} contains the unknow n degrees of freedom at the nodes of the elements. Equation (12) 
is solved using an asymmetrical, parallel, direct sparse solver to obtain the unknow n nodal parameters, { d>} . The solver 
uses a compressed column storage scheme to reduce storage overhead* Only the nonzero coefficients within [A] are 
stored, along with two pointer arrays that store the column numbers and starting indices of these nonzero coefficients* 
To achieve efficient solutions, the sparse solver uses two accelerators : equation reordering to reduce fill during the 
factorization of [A] and parallelization (i.e., the equation solver runs on multiple processors simultaneously). The finite 
element solution vector gives an approximation for the acoustic wall pressure opposite the test liner, for a 

specified wall impedance function, twaij. To avoid approximating the measured lower wall acoustic pressure, [zi ) , 
on the finite clement grid, the finite element grid is chosen so that each measurement point, z = Zh corresponds to a 
grid line in the finite dement analysis. The finite element approximation to the acoustic wall pressure opposite the 
lest liner, PrimU/), dong with the measured values, Pm(zi)* provides the necessary information for the impedance 
eduction technique discussed in the following section. 


VI, Impedance Eduction Technique 

G iven an initial estimate, £ wa n = timdah f° r the wall impedance function, the finite element methodology is used 
to provide the values for the acoustic pressure [i.e., /?fem(£/) I d measurement points on the wall opposite the lest 
liner. Values of ^initial are then iteratively updated until a wall impedance function, t wa i u is found that reproduces the 
measured wall pressure to within a specified tolerance. It is necessary, however, to automate the iterative procedure so 
that each new estimate f or Cinitifii is an improvement over its predecessor* We therefore implement an automated search 
procedure using an optimization algorithm* The objective function, (|>(t W aii), for the optimization algoritlun is defined 
as the differences between the measured acoustic pressure on the wall opposite the test liner and that computed by the 
finite dement method: 7 

nwatl 

<K£wau) = ^ I pM(fr) “ £FEm(2/) I (13) 

where nwall is the number of measurement points on the wall opposite the test liner. The optimization algorithm 
returns the value of the impedance function that minimizes the above objective function. It is a simple matter to show 
that the global minimum of the wall objective function, ^(^ wa u), is the unknown impedance function of the test liner. 

One of the most important aspects of the impedance eduction technique is the optimization algorithm. The opti- 
mization algorithm chosen was Stewart’s adaptation of the Davidon -Fletcher Powell (SDFP) optimization algorithm 
that uses a finite difference approximation to the gradient of the objective function. 8 This gradient- based optimizer 
has the disadvantage that it may converge to local optima (if they exist) and may become stuck in Lhat portion of the 
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impedance space where the objective function is extremely flat. However, these shortcomings are more than mitigated 
by the fact that SDFP lends to converge faster and give more accurate results thEui many of its competitors. Further, 
the occurrence of multiple local optima has to date not been observed with the impedance eduction process. However, 
since we arc not aware of any definitive proof that multiple local optima cannot exist for this process, we continue Lo 
use multiple initial estimates CCtmtiat) to strengthen our confidence that the global optimum is achieved. 

Recently, we have added a globed -based optimizer, die “Genetic Algorithm to our suiLc of optimization tools. 
Although the Genetic algorithm 9 locates a global optimum, it tends to he much more computationally expensive 
because it performs a global search in the impedance plane. It was therefore used only sparingly in this stud)’, and 
then only to spot-check the integrity of impedances educed using SDFP in the vicinity of where the anomaly occurs. 
In all cases, no noticeable differences between the impedance educed using SDFP and the Genetic Algorithm were 
observed. Therefore, only the SDFP results tire presented, 

VII. Results and Discussion 

I N this section, the impedance eduction methodologies are tested first on data synthesized from an exact mode 
solution and afterwards on measured data. r lhese two examples illustrate whether the anomalies occur both w ith 
synthesized and measured data. Impedances are educed at twenty-seven excitation frequencies ranging from 0.4 
to 3.0 kHz in 04 kHz increments. The impedance eduction codes use F90 Fortran with double precision (i.e., 64 
bit) arithmetic Eind were designed to run in a multiprocessor environment using a shared memory implementation. 
This approach has been chosen because the primary hardware to be utilized was a Silicon Graphics, Altix 3700 
system. In the current implementation the resistance, 0, mid reactance, °f ^e finer arc the design variables 
Kwall — 0 waU 4- j'Xwaiih The optimizer, SDFP, runs only in scquenlkd mode, using central finite differences to compute 
the gradient of the objective function, ^(t wa u). The finite difference step size is I x 10 _K and a stopping criteria of 
<K?waJl) <1x10 R was used lo terminate the search. Hie optimizer was constrained to realistic resistance values, 
0 VV aii >0.0 and the duct geometry is that of the NASA Langley GIT 2 (i.e.. Fig. 1) for which L ~ 0.813 meters, L\ “ 
0.203 meters , and L 2 =0.610 meters. The GIT has a square cross-section, so that the width (i.e., the distance from 
sidewall to sidew all) equals the duct height (i.e,, the distance from bottom to top wall), llie height of the GIT is 0.051 
meters 


A. Synthesized Data 

I11 the first example, we use a liner 0,407 meters long (L 2 — L\ ~ 0.407 meters) without a faceslieet mass reactance. 
The specified impedance spectrum is 

^walt ~ ^ wall i COt ( kd ) ( 14) 

where a low resistance of 0. 15 (8 wa u = 0. 15) is used and the depth of the cavity, d, is chosen such that the resonance 
frequency of the liner is at 1 .5 kHz. These liner parameters were chosen because they produce an impedance spectrum 
similar to that of liners for which anomalies have occurred based on data acquired in the GIT. The input data used 
for this first example is synthesized from the exact mode solution for outgoing waves in an infinite duct using the 
liner model of Eq. (14). To simulate an infinite duel without reflections and a single outgoing mode, we remove the 
(nirdwaJL duel sections upstream mid downstream of the liner mid place the source mid exil planes at z = L\ mid " == Z#, 
respectively In the presence of the hard Eind soft wall sections upstream mid downstream of the liner, higher-order 
modes mid reflections tire generally present in the vicinity of near fields generated by the leading L\) mid trailing 
edge (" = L 2 ) of the test specimen. For this reason the hard wall sections have been removed. The three sets of data 
needed lo educe the impedance spectrum given in Eq. (14) using the CHE model Eire the source pressure {/?,}, the 
lower wall acoustic pressure [pm(zi))* and the exit impedance, £ e xh- These three sets of data were obtained from the 
mode solution 

00 

p(z,x) = y A,,cos(X„x)e~ ,fC " z 

n= 1 

where K n Eire complex zeroes of the transcendental function 

E = -k n lan(X„H) 

<= wall 

M m K„ \ + AV/avcA-V \ 2 



(15) 

(16) 
(17) 
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and the eigenvalue, K „ , and axial propagation constant, K„, are related by the dispersion relation 


Kn 

k 


Jl-tl-A&e) [A* A) 


( 18 ) 


Note that for a given value of X n * the axial propagation constant K, : satisf ying Yl q. (18) lias two roots. Here, we 
will choose the root that corresponds to right-moving waves in the duct. It is easily shown that right-moving waves 
are identified as those modes for which the axial propagation constant, K ny has a zero or negative imaginary part 
and a positive real part, respectively. The method of analysis described above requires a calculation of the mode 
eigenvalues, that are the complex zeroes of the transcendental function, E , The method used to obtain these 
roots w as the New ton-Raphson iterative scheme. From the mode solution in Eq. (15), the data required to educe the 
impedance spectrum in Bq. (14) was obtained. Only the least attenuated mode was allow ed to propagate through the 
duct and the mode coefficient, w as set al a 120 dB level with zero phase. A uniform spatial grid with 129 points 
in the axial direction (N = 129) and 21 points in the transverse direction (A/ = 21) was used to educe the impedance 
from the synthesized data. 

Table 1 

Average Mach number and mean static conditions for synthesized data 



Pa 

To 

0.000 

101425.45 

297.48 

0.085 

10 1473.94 

296.48 

0.253 

101453,30 

295.15 

0.401 

10121 1.70 

293.37 


The impedance spectra was educed (using the CUB model) for the four uniform flow Mach numbers, mean static 
pressures (in Pascals), and mean static temperatures (in Kelvin) given in table 1. The resistance mid reactance spectra 
educed using the CHB model are compared to the exact spectra [see Bq. (14)] in Figs. 2 mid 3 t respectively. To the 
level of accuracy presented in Figs, 2 and 3, the educed iniped;inec reproduce the exact impedance in Bq. (14) without 
any anomalies. We have also noticed that the small values of the finite difference step sizes and stopping criterion 
used in SDFP were also important parameters contributing to such excellent agreement betw een the educed and exact 
spectrum, 

A major drawback introduced by the local, gradient-based optimizer (i.e., SDFP) is that it may converge to a local 
minimum if local minima exist. In this situation, different design points (i.e., educed impedances) may be obtained 
if the optimizer is initialized from different points in the design space. Because of reduced turnaround time resulting 
from usage of the parallel, sparse solver, an exhaustive evaluation of the effects of “starting location” (i.e., initial 
impedance) could be conducted. Several distinct starting locations were used ;ind no local minima were found. 

R. Measured Data 

In the second example, input data required for the impedance eductions is obtained from measurements taken in the 
NASA Langley GIT. 2 As described earlier, the CHE and LEE are 2D impedance eduction models. Thus, only data 
that is acquired at excitation frequencies below the “cut- on” of higher- order modes in the hard- wall section of the 
GIT should he expected to he properly described by these models. Due to the sound absorbing properties of the Lest 
liner, the possibility of higher-order modes cannot be avoided in the lined region. In Lhc presence of the liner, higher- 
order modes and reflections arc generally present in the vicinity t)f near fields generated by the leading (z = Li) £ind 
trailing edge {z — Li) of the test specimen. The test liner is a conventional pcif orate type liner composed of a punched 
dunlin um face sheet bonded to a honeycomb core that was in turn bonded Lo a rigid back plate. The face sheet is 
0,635 x 10' 3 meter thick with holes that w ere 0.991 x 10" 3 meters in diameter. The open area ratio is 0,087. The 
honeycomb core has a cell size of 9,525 x 10 -3 meters and a depth of 38400 x 10“ 3 meters. Because the test liner 
w r as designed from homogeneous material, the impedance was assumed to be invariant over the length of the liner. 
The incident sound pressure level (SPL) is kept at a low 7 excitation level (approximately 120 dB) to minimize liner 
nonlineari ty effect s . Finally, the decision w as made to present only zero flow 7 results using the measured data. This is 
done to avoid processing measured data (obtained in the presence of a boundary layer) through impedance eduction 
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codes that neglected the effects of the boundary layer. Flow results in the presence of measured data are currently 
being analyzed with the LEE model and will be reported at a later dale. 

Figure 4 compares the educed normalized impedance obtained from the CHE model using the baseline grid (i,e., 
N — 129, mid M =21) with that measured in the I jmglcy Rcsciirch Center Normal Incidence Tube (NJT). The educed 
normalized rcsisOmec and reactance spectra sire plotted using a dual axis system, with the normalized resistance mid 
reactance referenced lo ihe left and right axis, respectively. The NTT results sire used as the baseline against which to 
measure the accuracy of the GIT results because they are normally considered more accurate in the absence of grazing 
flow. However* NIT measurements cannot be made in the presence of g razing flow so GIT results are often relied 
upon. As expected, both the NIT and GIT normalized reactance show a cot (fa/) type of behavior and the normalized 
resistance has a low value of approximately 0.15 (see Fig. 4). The NIT and GIT spectra are in very good agreement 
except at 0.4, 0.5, 2.8, and 3.0 kHz. For example, the educed normalized resistance at 0,4 kHz is lower than that 
measured in the NIT and the value at 0.5 kHz is high (i.e., when compared to NIT measurements). Subsequently, 
the slope of the edueed resistance curve in the region 0.4 to 0.5 kHz is opposite in sign to that measured in the NIT. 
Further, the change in slope of the educed normalized resistance spectrum near 2,8 kHz mid rapid rise in its slope 
near 3.0 kHz is not only surprising, but is not reproduced by the NIT measurements, This strange behavior at low 
frequency (0.4 and 0.5 kHz) mid high frequency (2.8 mid 3.0 kHz) has been observed in numerous other tests with 
different test liners mid is referred to here as impedance eduction anomalies. Note that the anomalous behavioral 0,4 
and 0,5 kHz (see Fig. 4) occurs in the vicinity of the anti -resonant frequency on the Low frequency end of the spectra, 
and the anomalous behavior at 2.8 and 3.0 kHz occur in the vicinity of the anti -resonant frequency on the high end of 
the spectra. 

To determine if ihe anomalous behavior observed iu Fig, 4 was due to lack of grid refinement or a breakdown of 
the CUE model, the impedance spectrum was educed using the LEE model and a grid refinement study was performed 
on both models. Figure 5a shows the results of three different grid resolutions using the CHE model. Again the 
resulLs fire show n on a dual axis system w ith the normalized resistance mid reactance reference to the left mid rigliL 
axis, respectively. The baseline grid, GO, has 129 evenly spaced points (N = 129) along the axis of the ducL mid 21 
(Af =21) evenly spaced points in the vertical direction (i.e., betw een the top mid bottom wall). The density of the 
second grid, G2, is obtained by doubling the density of the baseline grid. Finally, the third grid, G4, is obtained by 
doubling the density of the G2 grid. Thus, the G4 grid is four times denser (in each direction) than the baseline grid. 
Figure 5b shows results of the grid refinement study on the LEE model* Both models are observed to have converged 
on the baseline grid for this impedance spectrum* 

Figure 6 compares (using the dual axis system) the educed normalized impedances on the baseline grid for the 
CUE (Fig, 5a) and the LEE (Fig, 5b) models* Both models give educed impedances that are nearly identical (at least 
on the scale for which the results are plotted) and both models produce the anomalies at 0,4, 0.5, 2.8, and 3,0 kHz, 
When one considers that 

1. each model supports a different set of equations, 

2. each model uses a different approximation to the acoustic pressure, and 

3. each set of results are obtained using the SDFP optimization algorithm* 

then the excellent comparisons betw een the models in Fig. 6 are indeed gratifying. Whereas the comparison between 
CHE and LEE in Fig. 6 is not mi infallible basis for judging the accuracy of these models, it does show that the 
impedance eduction theory (as applied here) is not model dependent. Such a result is always a necessary step in 
assessing the source of the impedance eduction anomalies* 

One potential explanation of the enigmatic behavior at frequencies close to an anti -resonant frequency (see Fig. 4) 
is the failure of ihe SDFP optimization algorithms. Figure 7a shows a contour map of the wall objective function, 
at l h e lowest frequency (i.e., 0,4 kHz) where the anomaly occurs. The contour map was constructed using 
Ihe CHE model and on the baseline spatial grid, GO. Here, 81 evenly spaced points along the normalized resistance 
axis and 81 evenly space points dong the normalized reactance axis were used to construct the contour. The center 
of the eye of the contour (i.e,, the optimum point) is located at twaii = 0.07 — 3.72/. This corresponds to the value of 
Ihe optimum impedance obtained using SOFT (Fig. 4). The diameter of the eye is approximately 0,1 dimensionless 
units mid is larger than the edueed normalized resistance of the test liner. Because the edueed normalized impedance 
is the location of the center of the eye, then for the current low resistance liner, the unknown normalized impedance 
cannot he determined with belter accuracy using the measured data. This conclusion is also supported by Fig, 7b Lhat 
compares the measured lower wall SPL (in decibels) and phase (in degrees) w ith that predicted by the CHE model. 
The results are plotted on a dual axis system with the lower wall SPL mid phase reference to the left and right axis, 
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respectively. The symbols in the figure tire located at the 31 lower wall measurement points (i.e., z =■ Z/). Recall that at 
the optimum point, the lower wall SPL and phase returned by the CHE model should match those measured in the GIT. 
As shown in Fig, 7b, the lower wall SPL and phase predicted from the CHE model are well matched Lo the measured 
v ill ucs. Although, Lhc results tire not shown for the sake of brevity, similar trends to that in Fig, 7a and Fig. 7b were 
observed at 0.5 kHz. 

As illustrated in Fig. 4, the discrepancy between NIT and GIT impedances are more prominent at the two highest 
frequencies (2.8 and 3.0 kHz). Figure 8a shows a contour map of the wail objective functions at 2,8 kHz. The eye of 
the contour is more diffused than at 0.4 kHz and the diameter of the eye has doubled from that observed at 0,4kIIz (see 
Fig. 7a). The optimum impedance from the contour map matches that from SDFP (see Fig. 4) quite well. Figure 8b 
shows a comparison of the measured lower wall SPL and phase (using the dual axis system) to that predicted from 
the CHE model. Although, the measured phases are well matched, significant differences are observed between the 
measured SPL and those predicted by the CHE model. These SPL differences are more prominent upstream of the 
leading edge of the liner (i.e., z < 0.203 meters) mid in the vicinity of the near field produced by the leading edge of 
the liner. 

Figure 9a shows a contour map of the wall objective function at 3.0 kHz. Differences between the normalized 
impedance measured in the NIT and that educed using the CHE model are the most prominent at this frequency. The 
eye is strongly diffused with a diameter of nearly 0.6 dimensionless units. A comparison of the measured lower wall 
SPL and phase with that predicted by the CHE model is given in Fig, 9b, The differences between the measured and 
predicted SPLs mid phases are the most noticeable upstream of the leading edge of the liner (i.e., z < 0.203 meters) 
and in the vicinity of the near field produced by the leading edge of the liner. Much smaller eyes were observed at 
frequencies away from an anti -resonant frequency (i,e., where the anomalies do not occur). For example. Fig. 10a 
shows a contour map of the wall objective function at 2.5 kHz, This contour map is representative of what was 
observed at frequencies where the anomalies did not occur. Note that, in this ease, the eye of Lhc contour appears 
lightly wound with a much smaller diameter (approximately 0.03 dimensionless units) than is observed at frequencies 
close Lo mi anti-resonance (i.e., 0.4, 0.5, 2.8, and 3.0 kHz). Figure 10b compares the measured lower wall SPLs mid 
phases with that produced by the CHE model. Note that the measured SPLs and phases arc well matched to those 
produced by the CHE model. 

A brief study was conducted to investigate the potential sensitivity of the educed impedance to errors in the exit 
impedance. Of the input data parameters, the exit impedance boundary condition is of concern, at least in part because 
it is determined from acoustic pressure measurements with microphones not included in the objective function [see 
Eq, (13)]. Figure 1 1 show the educed normalized resistance obtained from the CIIE model by decreasing the normal- 
ized exit resistance by 0.1 over the measured value (the curve with triangle symbols) mid by increasing the measured 
normalized resistance by 0. 1 (the curve w ith circles). TTiis perturbation value was chosen because resistance variability 
of this amount has been observed in practice at selected frequencies. Note that at mid around the frequencies close to 
an anti -resonant frequency (i.e., 0.4, 0.5, 2.8, mid 3.0 kHz) measurable changes in educed resistance are observed. 

Having identified errors in the exit inipediinec as one potential source of the anomalies, two remedies arc suggested 
to minimize this error source: 

1. Allow the optimization algorithm to educe the normalized exit impedance by including Lhc normalized exit 
resistance and normalized exit reactance as design variables in the optimization algorithm. 

2. Replace the exit imped mice boundary condition with an acoustic pressure boundary condition [Eq. (5)], such 
that the data used in the boundary condition are not independent from the remainder or the data used in the 
impedance eduction process. 

Figure 12a shows the educed normalized resistance When each of the above remedies is implemented in Lhc GHE 
model. Also plotted in the figure tire the normalized resistances measured in the NIT mid those educed using the 
exit impedance measured in the GIT. Note that the anomalies at 0.4 and 2.8 kHz are reduced substantially (i.e., the 
educed resistances are much close] - to the NIT measurements) when either of the suggested remedies are implemented. 
At 3.0 kHz, the suggested remedy of educing the exit impedance gives results closer to the NIT measurements than 
implementing the exit pressure boundary condition. On the other hand, at 0,5 kHz, implementing the exit pressure 
condition brings the educed resistance closer to the NIT measurements than the educed exit impedance approach, 
although, the improvement is quite modest. Figure 12b shows the results for the educed normalized reactance. There 
is very little effect of either remedy on the educed reactance except at the highest frequency (3,0 kHz), At 3.0 kHz, 
implementing the exit pressure condition reduces the high value of the educed normalized reactance obtained by using 
the measured exit impedance by 0.35 dimensionless units, whereas the reduction is only about 0,15 units using the 
educed exit impedance. Generally speaki ng, implementing the exit pressure boundary condition is observed to reduce 
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the normalized impedance anomalies more than educing the exit impedance (Fig, 12a and Fig. 12b). Figure 13a shows 
the contour map when the exit pressure boundary condition was implemented. Note that the diameter of the eye is 
now only fifty percent of the value obtained when the measured exit impedance was used Lo construct Lhe contour map 
(see Fig. 9a). Figure 13b shows a comparison of the measured lower wall SPLs and phases to that predicted by the 
CHE model when the measured exit impedance is used mid when the exit pies sure is used as the terminating boundary 
condition. Note that there is improvement in the SPL and phase comparisons downstream of the trailing edge (he,, 
Z > 0.610 meters) of the liner but little improvement upstream of the leading edge (i.e., z < 0,203 meters). 

While the impedance educed at 3.0 kllz is improved by the exit pressure boundary condition, the educed impedance 
at tills frequency remains an outlier in the educed impedance spectrum (i.e., w hen compared to the NIT measurements). 
Observe that 3.0 kllz is very near anti “resonance, and is also approaching cut-on of the next higher-order mode. Since 
the liner attenuation is quite minim al near anti -resonance, differences in impedance such as that observed between 
the NTT and GIT results in Fig. 4 has traditionally not been of great concern. Further, it should ill so be noted that it 
near-evanescent mode is generated near the leading edge of the liner at a frequency approaching cut-on (e.g., 3,0 kHz), 
and this mode may survive to the source plane of the computational domain, because of it slow axial decay and the 
absence of acoustic treatment upstream of the leading edge of the test liner. This would result in a source plane that is 
not planar in the transverse direction (i.e., the x direction) as assumed in the CHE and LEE models. The propagation 
of a high -order, uear-evancscent mode, upstream to the source plane is one plausible explanation for the discrepancy 
between measurement and theory at the frequency closest to cut-on. One possible remedy for tins would be to use the 
acoustic pressure profile across the source plane as the source boundary condition in the impedance eduction models. 
The GIT’ test window is designed to acquire these data, and this option may be pursued in future endeavors, 

V1IL Conclusions 

ased upon the results of this study, the following preliminary conclusions can lie drawn: 

1. In the current studs - , anomalous impedances were educed only at frequencies near anti -resonance or cut-on of 
higher-order modes. 

2. Sensitivity of the educed impedance to small variations in the measured exit resistance is increased for these 
frequencies where anomalous imped mice s are educed. 

3. This anomalous behavior occurs at identical frequencies for cadi of the impedance eduction models, whether 
based on Lhcconvecled Helmholtz equation or the linearized Filler equations. 

(a) Grid refinement studies indicate the choice of grids is not the cause of these anomalies. 

(b) Two optimization algorithms - one local and one global - provide similar anomalous results. 

4. These anomalies are reduced by replacing the exit impedance boundary condition with an acoustic pressure 
boundary condition at the exit plane. 

The results of tliis study offer evidence that, at least for the no-flow condition, it is important to properly model the 
termination boundary condition for the computation domain. It appears particularly important to note that the common 
usage of an an echoic, plane wave termination boundary condition in ducts where the “true” termination ts reflective, 
may act as a trigger for anomalous results. The effects of mean flow on this “conclusion” need to be evaluated and tins 
effort is currently underway. 
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I'ig 2. Exact anil educed resistance from CHE model using synthesized data: shows no anomalies. 
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Fig 3. Exact and educed reactance from CHE model using synthesized data: 


shows no anomalies. 
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Frequency, kHz 


Fig 4. Comparison of the educed normalized impedance from the CHE model using measured data to normal incidence 
lube (NIT) measurements without flow: shows anomalies at frequencies close to anti -resonance {0,4, 0,5, 2.8, and 3,0 
kHz). 
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Fig 5a. Educed normalized impedance spectra produced by Lhe CHE mode] on llircc grids (GO, G2, and G4) using 
measured data and no flow: shows the invariance of the educed normalized impedance spectra to the change in grid 
resolution. 


14 of 29 



Frequency, kHz 


Fig 5b, Educed normalized impedance spectra produced by the LEE model on three grids {GO, G2. and G4) using 
measured data and no flow: shows the invariance of the educed normalized impedance spectra to change in grid 
resolution. 
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Fig 6. Comparison of educed normalized impedance spectra obtained from the CHE and LEE model using measured 
data on the baseline grid, GO, without flow: shows convergence of the CHE and LEE models. 
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Fig 7a, Contour map of the wail objectiv e function produced by tlie CI1H model at 0,4 kHz using measured data and 
no flow: shows that the optimum point of the contour map equals that produced by SDFR 
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Fig lh. Comparison of the measured lower wall sound pressure level (SPL) and phase at (J.4 kHz to that predicted by 
the OIL model using measured data and no How; shows a good match of measured and predicted SPL and phase. 
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Fig Sa. Contour map of the wail objective function predicted by the CIIH model at 2.8 kHz using measured data and 
no flow: shows a moderate size eye. 


19 of 29 



130 i- 


► 

o 


128 - 


Measured SPL 
Predicted SPL (CHE model) 
Measured phase 
Predicted phase (CHE model) 


- 180 



- 90 


- 0 


- -90 


0.2 0.4 0.6 

Axial coordinate z, meters 


-180 


Fig 8b. Comparison of the measured lower wall sound pressure level (SPL) and phase at 2.8 kHz to dial predicted by 
die CHF model using measured data and no Bow: shows SPL errors upstream of die leading edge of the liner and in 
die vicinity of the near field produeed by the leading edge of the liner. 
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Fig 9a. Contour map of the wall objectiv e function predicted by the CHE model at 3.0 kHz using measured data and 
no flow: shows an extremely huge eye. 
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Fig 9b. Comparison of the measured lower wall sound pressure level (SPL) and phase at 3,0 kHz to that predicted by 
the CHB model using measured data and no flow: shows SPL and phase errors upstream of the leading edge of the 
liner mid in the vicinity of the near field produced by the leading edge of Lhe liner. 
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Fig 10a. Contour map of the wall objective function predicted by the CUE model at 2.5 kHz using measured data and 
no flow: shows an extremely small eye. 
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Fig 10b. Comparison of llie measured lower wall sound pressure level (SPL) and phase at 3.0 kHz to that predicted by 
the CHE model using measured data and no flow: shows a good comparison of SPL and phase. 
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Fig 1 L Effects of exit impedance resistance errors of ±0.1 on the educed resistance using the CHE model, measured 
data, and no flow: shows sensitivity of the educed normalized resistance to measured exit resistance errors aL and near 
the anomalous frequencies 0.4, 0.5, 2.8, mid 3.0 kHz. 
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Fig 12a. Educed normalized resistance using ihc exit pressure, the exit impedance, and the educed exit impedance 
boundary condition using the CHE model, measured data, and no flow: shows a reduction in the anomalies aL 0.4, 2.8, 
and 3.0 kHz using the educed impedance mid exit pressure boundary condition. 
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Fig 12K Educed nonnalized reactance using the exit pressure, the exit impedance, and the educed exit impedance 
boundary condition using the CHE model, measured data, and no flow: shows a reduction in the anomaly at 3.0 kHz 
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Fig 13a. Contour map of the wall objective function predicted by the CHE model at 3.0 kHz. using measured data, the 
exit pressure boundary condition, mid no flow: shows a fifty percent reduction in the size of the eye compared to the 
exit impedance boundary condition (see Fig. 9a). 


28 of 29 


Measured SPL 

Predicted SPL (CHE model) exit impedance 
Predicted SPL (CHE model) exit pressure 
Measured phase 

Pedicted phase (CHE model) exit impedance 
.^Predicted phase (CHE model) exit pressure 


130 

r o 


> 

128 

_ < 


H ■ ■ ■ ■ ■ 

126 


~ 124 
CD 

■o 

£L 122 

CO 1 


120 


118 




116 



- 270 


- 1 80 


- 90 



- -90 


J =1 -180 


0.2 0.4 0.6 

Axial coordinate z, meters 


0.8 


Fig 1 3b. Comparison of the measured lower wall sound pressure level (SPL) mid phase at 3.0 kHz with that produced 
by the CHE model using measured data, the measured exit impedance, the exit pressure boundary condition, mid no 
flow : shows better comparison to measured data using the exit pressure boundary condition. 
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